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Derivation d'un modele de type Saint-Venant 
non hydrostatique; 
Comparaison avec les modeles de type 
Boussinesq 

Resume : A partir dcs equations de Navier-Stokes a surface libre, on obtient 
deux modeles moyennes sur la verticale, non hydrostatiques qui etendent le 
systeme de Saint-Venant et incluent le frottement et la viscosite. La complexite 
des formulations obtenues depend du niveau d'approximation retenu pour la 
pression du fluide. Les modeles obtenus sont compares aux formulations de 
type Boussinesq. 

Mots-cles : Equations de Navier-Stokes, equations de Saint-Venant, equations 
de Boussinesq, surface libre, termes dispcrsifs 
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1 Introduction 

Despite the available numerical results obtained by the simulation of the Navier- 
Stokes equations, there exists a demand for models of reduced complexity such 
as shallow waters type models. 

Non-linear shallow water equations model the dynamics of a shallow, rotat- 
ing layer of homogeneous incompressible fluid and are typically used to describe 
vertically averaged flows in two or three dimensional domains, in terms of hor- 
izontal velocity and depth variation, see Fig. [TJ This set of equations is par- 
ticularly well-suited for the study and numerical simulations of a large class of 
geophysical phenomena, such as rivers, coastal domains, oceans, or even run-off 
or avalanches when modified with adapted source terms [7J. 

The classical Saint- Venant system [3] with viscosity and friction (T4l [T71 [T3] 
is well suited for modeling of dam breaks or hydraulic jump but due to the 
hydrostatic assumption it is not well adapted for the modeling of gravity waves 
propagation. 

For the modeling of long wavelength, small amplitude, gravity waves, the 
Boussincsq system [§1 EH] is used. The Boussinesq equations are obtained 
from the Euler equations i.e. ignoring rotational and dissipative effects [H lilt 
H2 HH [201 123]- In practice, the use of such models ignoring rotational and 
friction effects at the bottom may be very restrictive. Furthermore, even when 
well posed, the Boussinesq models often exhibit a lack of conservation energy 
that is odd since they are derived from Euler equations [SI [5] • 

The objective of this paper is twofold. First, we want to extend the Saint- 
Venant system so that the long waves propagation can be modeled and second 
we aim at comparing/unifying the obtained formulation with the Boussinesq 
system, see Fig. [U The paper is organized as follows. In section [2j we recall 
the Navier-Stokes system with a free moving boundary and its closure. We also 
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Navier-Stokes equations 
for incompressible free surface flows 



Shallow water assumption 



^ Hydrostatic 



models 



Saint- Venant system 



Non hydrostatic models 



Boussinesq system 



Multilayer Saint- Venant system Extensions of the Saint- Venant system 

Figure 1: Averaged models derived from Navier-Stokes equations. 



present the Saint- Venant and Boussinesq assumptions and the associated rescal- 
ing. In section [3] we recall the Shallow Water system and show the hydrostatic 
Boussinesq system assumption corresponds to the classical Saint- Venant system. 
In section^ the hydrostatic assumption is relaxed and we obtain two formula- 
tions of growing complexity extending the Saint- Venant system and depending 
on the level of approximation chosen for the fluid pressure. 



2 The Navier-Stokes system 

Let start with the Navier-Stokes system [TB] restricted to two dimensions with 
gravity in which the z axis represents the vertical direction. For simplicity, the 
viscosity will be kept constant throughout the paper. Therefore we have the 
following general formulation expression: 

du dw 

du | ^du [ ^du ^ 1 dp _ dT, xx dT, xz 

dt dx dz p dx dx dz ' 

dw dw dw 1 dp dY, zx dY, zz 

dt dx dz p dz dx dz 

and we consider this system for 

t > to, x £ R, Zb(x,t) < z < rj(x,t), 

where r](x,t) represents the free surface elevation, u = (u,w) T the horizontal 
and vertical velocities. The water height is H = r\ — Zb, sec Fig. [2] We consider 
the bathymetry Zb can vary with respect to abscissa x and also with respect to 
time t. The chosen form of the viscosity tensor is 

2-, xx = 2v—, Y, xz =v{— + —), 
dx v dz dx' 

dw ,du dw, 

2- lzz = 2v—, E zx = i/f— + —J, 
dz dz dx 
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with v the viscosity coefficient. For a more complex form of the viscosity tensor 
using eddy and bulk viscosities, the reader can refer to [15j . 




Free surface 



H(x,t) 



<pb( x tt) 



Bottom 



Figure 2: Notations: water height H(x,t), free surface r](x,t) and bottom 
Zb(x,t). 



2.1 Boundary conditions 

The system ([I])-© is complete with boundary conditions. The outward and 
upward unit normals to the free surface n s and to the bottom xib are given by 



'i+(ir v 1 7 v i+ (^ 

Let T,t be the total stress tensor with 

S r = --pl d 4 
P 

2.1.1 At the free surface 

Classically at the free surface we have the kinematic boundary condition 

- +Us -- Ws = 0, (4) 

where the subscript s denotes the value of the considered quantity at the free 
surface. Considering the air viscosity is negligible, the continuity of stresses at 
the free boundary imposes 

p a 

S T n s = n s , (5) 

P 

where p a = p a (x,t) is a given function corresponding to the atmospheric pres- 
sure. Relation ^ is equivalent to 

p a 

n s .E T n s = , and t s .S T n s = 0, 

P 

t s being orthogonal to n s . 
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2.1.2 At the bottom 

Since we consider the bottom can vary with respect to time t, the kinematic 
boundary condition is 

dz b dz b 

= (6) 
where the subscript b denotes the value of the considered quantity at the bottom 
and (x, t) i — ► Zb(x, t) is a given function. Note that Eq. © reduces to a classical 
no-penetration condition when z b does not depend on time t. 

For the stresses at the bottom wc consider a wall law under the form 



S T n fc - (n b .E T n b )n h = K(v b ,H)v b , (7) 

with Vf, = u b — (0, ^§f) T the relative velocity between the water and the bottom. 
If k(vi,,H) is constant then we recover a Navier friction condition as in [14j . 
Introducing laminar fc; and turbulent k t friction, wc use the expression 

«(vb,-ff, v) = h + k t H\v h \, 

corresponding to the boundary condition used in |17j . Another form of K(vb, H) 
is used in [7J and for other wall laws, the reader can also refer to [IB]. Due to 
thermomechanical considerations, in the sequel we suppose K(vb,-ff) > and 
K(vb, H) is often simply denoted k. 

Let t(, satisfying t^.n^ = then when multiplied by tf, and n b , Eq. ([7]) leads 

to 

t^.E^nh = Kv b .t b , and v b .n b = 0. 



Remark 1 If the boundary condition |[7J) was written under the form YiT-n b = 
K(vb,-ff)vf, as in Ferrari et al. \ltft Eq. (2.25), p. 217], then in absence of 
friction and viscosity, this would give 

Pb = 0, 

that is not correct. 

2.2 The rescaled system 

The physical system is rescaled using the quantities 

• h and A, two characteristic dimensions along the z and x axis respectively, 

• a s the typical wave amplitude, a b the typical bathymetry variation, 

• C = \fgh the typical horizontal wave speed. 

Classically for the derivation of the Saint- Venant system, we introduce the small 
parameter 

h 

£ =\- 

When considering long waves propagation, another important parameter needs 
be considered, namely 
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and wc consider for the bathymetry 5£ = O(S). Depending on the application, 
(5 can be considered or not as a small parameter. For finite amplitude wave 
theory and assuming Zb(x,t) = z®, one considers e -C 1, 6 = 0(1) whereas the 
Boussinesq waves theory requires 

e«1 and U r = 0(l). 

where U r is the Ursell number defined by U r = p-, see [33]. All along this 
work, we consider eCl whereas, even if the parameter 5 is introduced in the 
rescaling, the assumption S <C 1 is not considered (paragraphs 14. 1[ 1 4.21 and |4~3]) 
except when explictly mentioned. 

As for the Saint- Venant system [HI \T7\ , we introduce some characteristic 
quantities : T = X/C for the time, W = a s /T = a^/T = eSC for the vertical 
velocity, U = W/e = SC, for the horizontal velocity, P = pC 2 for the pressure. 
This leads to the following dimensionless quantities 

-X_Z_T]~t 

A ' K 1 ' a s ' T ' 

p _ u w 
p = —, u = — , and to = — . 
y P' U' W 

Note that the definition of the charateristic velocites implies 5 = ^ so S also 
corresponds to the Froude number. When S = 0(1) we have U ~ C and 
we recover the classical rescaling used for the Saint- Venant system. For the 
bathymetry z\, we write Zb(x, t) = Z^(x) + b(t) and we introduce Zb = Zb/h and 
b = b/ab- This leads to 

dz b db db dz b dz b 

— — = ecu— = W—, and — — = e— — . 
at at at ox ox 

The different rescaling applied to the time and space derivatives of Zb means that 
a classical shallow water assumption is made concerning the space variations 
of the bottom profile whereas we assume the time variations of Zb lie in the 
framework of the Boussinesq assumption and are consistent with the rescaling 
applied to the velocity w. 

We also introduce v = and we set K= ^. Note that the definitions for 
the dimensionless quantities are consistent with the one used for the Boussinesq 
system [50] [H] . Note also that the rescaling used by Nwogu [TH] differs from 

2 

the preceding one since Nwogu uses w = j^w. 

As in [14|. I17j , we suppose we are in the following asymptotic regime 

v = sisq, and k = ekq, 

with kq — Ki t o + £Kt,o(v&, H) 7 kj.o being constant. 
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This non-dimcnsionalization of the system (d])-© leads to 
0, 



du dw 
dx dz 



c du .,.95 _ du dp 9 d / du 

£ S + e5 2 u— + eS 2 w— + e-£ = e 2 6— 2is — 
dt dx dz dx dx V dx 



d_ 

dz 



dii 2 
ov — +e bv 

dz 



e 2 S 



dw _ dw _ d r i 
dt dx dz 



d_ 

dx 



du 
ITT 

dz 



dw\ 






~di) 






9+ 


dp 
1 


= -1 


du> 
dx 


) + 


dz 



2v Q 



di 



(8) 



(9) 



(10) 



with the boundary conditions (JU), (O, ([!]) and ((TJ) becoming 



dfj dfj „ 
at ax 

2edv -z-r -p s 

dz „ 



2 2 9?7 ^ du 



ax V dz 



dw 
dx 



5v 



6vq e 



du 
dz s 

db _ dzb 

d! + Ub ^ 

dw 
dx 



2 dw 
dx 

w b = 

du 

+ 7F 
az 



5?7 / <9u 
ax \ ax 



dz b ( du 
dx \ dx 



-5p a 



dx 



Pb 



dz b ( dw 
dx \ dz 



dSb ( z du 
Pb - zvq— b — 
dx \ dz 



- . idhr- df, 

Ub+£ -M^-m) 



s 2 6 



dx 



(11) 
(12) 

(13) 

(14) 



(15) 



For the sake of clarity, in the sequel we drop the symborand we denote 



dt ■ 



3 The Shallow Water system 

In this section we first derive the expression of the fluid pressure p in the context 
of the Shallow Water assumption and then show the combination of the Boussi- 
nesq and hydrostatic assumption leads to the classical Saint- Venant system. 

The process used hereafter is similar to the technique employed by Gcrbeau 
and Perthame [14j to derive a formulation for the viscous Saint- Venant system. 
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3.1 The vertically averaged system 

Using the divergence free condition, the system ([8|)- (fT0|) is rewritten under the 
form 

du dw „ , „. 

^ + ^ = < 16 > 

c du p9 du 2 r9 duw dp , . d / <9u 
eS— + eS 2 —- +sS 2 — V e^- = e 8— 2v 



dt dx dz dx dx \ dx 

9 f c du 2 dw\ . 
8v —+e 2 8v —\, (17) 



dz \ dz dx 

dw duw dw 2 \ dp ^ 

dt dx dz J dz 

d ( c du - dw\ d ( dw\ 

+ d- x { £ ^d-z +£S ^dx-) +£S dz r°^) 

Due to the applied rescaling some terms of the viscosity tensor e.g. 

, .. <9 / dw 
£ 3 S— v 



dx \ dx 

are very small and could be neglected. But, as mentioned in [51 Remarks 1 
and 2] , the approximation of the viscous terms have to preserve the dissipation 
energy that is an essential property of the Navier-Stokes and averaged Navier- 
Stokes equations. Since we privilege this stability requirement and in order 
to keep a symmetric form of the viscosity tensor, we consider in the sequel a 
modified version of (|16 p -(|18 p under the form 

du ^ dw g 
dx dz 

r du r9 du 2 c9 duw dp 9 .9 / du 
£8— + eS 2 —- + e<r — h s^f = £ 8— 2v 



dt dx dz dx dx \ dx 



dz \ d 

9 c / dw rduw c .dw 2 \ dp d ( . du 

£ 8 — + 8—- + 8—- ) + = 1 -\- — — [ e8v 



d ( ' du s . 

'8v Q —\, (20) 



\ dt dx dz J dz dx \ dz 



corresponding to a viscosity tensor of the form 



„ 3« du dw 

^v—, Y, xz = T, zx = i/—, T, zz = 2v—. 
ox oz dz 



Remark 2 If we strictly follow Audusse JH Lemma 2.1], the chosen form of 
the viscosity tensor will not allow us to include under the form of a square term 
in the energy equality the quantity 

du dw 
dz dx 

But we will see in paragraph \4-S\ that due to the shallow water assumption, this 
quantity appear as a friction term. 
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From Eqs. flU}, it 

comes 



p s = Sp a + 2eS 



dw 

dz 



+ 0(e 2 S 2 ), 



so using Eqs. (|T3j) and (jT5J) one obtains 

az 



du 

c) z 



0(e), 



and an integration of Eq. (|2ip from Srj to z gives 

p - 6p a = 6rj- z + 0(eS), 

leading to 



(22) 



(23) 



dp 

dx 



O(S). 



The preceding relation inserted in ([20]) leads to 



and Eqs. 



and 



mean that 

u{x, z, t) = uq(x, i) + 0{e), 



(24) 



i.e. we recognize the so-called "motion by slices" of the usual Saint- Venant 
system. Then we introduce the averaged quantities 



u dz, u 2 = 

z„ <"7 - z b 



6ij 



i 2 dz, 



Si] - z b 

and the previous definitions involve 

u(x, z,t) = u + O(e), and u 2 = u 2 (x, z, t) + 0(e). 



(25) 



Note that the velocity u is exactly the one arising in the conservation law for 
the water height since an integration of Eq. (|19|) from z b to Srj with boundary 
conditions (|11[) and (Q3J) leads to 



dr) dz b d , s 



(26) 



with H$ = Srj — z b . Conversely an integration of Eq. (] 19|) from z b to z with 
boundary conditions (|11|) and (fT4")) leads to 



dz b d 
w = — r / Mdz 
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We use the approximations obtained in this paragraph to simplify the boundary 
conditions (|1 1|) -()15| ) and retaining only the high order terms we obtain 



dv . , dri 
_ +6Us -- Ws = 0, 



p s = 5p a + 2e5v 



dw 



8v Q 



du 

dz 



eS 



drj 
dx 



0(s 3 S), 

du 
dx 



Pi 



ox 



dz b dz b 



Sv 



_ dz b 
dx 



du 
dx 



Pb 



— E 



dz b 
dx \ 



2e5vq 



dw 



Pb 



^ dz b du 
dx dz 



3e 2 (dz^ 2 



s6 Ko [l + —[ — 



(28) 
(29) 

(30) 

(31) 



u b + 0{e 4 S). (32) 



Using the Leibniz rule i.e. 







dx 



b(x) 



aZ / 9 



a(x) 



f b[x) dg , db . . da 
dXl = L x) dx^^ 1 + d-x 9{a{x)) ~ d-x g{h{x)) > 



and the kinematic boundary conditions (|28p and (j3"Tj) , an integration of Eq. (|2TJ)) 
from z b to Sr/ shows that a solution to (fH?|) - (f2"Tj) satisfies 



^ J 9 r 1,2 2 



q nSr, 



Su 



and using Eqs. (f3"0"|) and (|32[) . we obtain 
d [ 5r > 



u + p) dz = e~5— 
dx 

2e<5^o 



d f 5r > „ d 



du 


-eS— 

g ' dx 


dz 


du 


+ e dZb 


~d~z 


b dx 



Zb 

du 
dx 
du 
dx 



2v —dz 
dx 



Pb 



d f 8n 
d dl 



n j 9 r i,2 2 \ , ,d r s \ d Uj 

I u dz + — (o u + p) dz — eo— Zv —dz 
J Zh vx J _ dx J z dx 



ax aa; 



2e<5^o 



9w 



cte dz 



3e 2 /SV 2 



*«o I 1 + — ( -j£ ) \ub + 0(e*8), 



(33) 



An expression for the pressure p can be obtained as follows. An integration of 
Eq. (|2"Tj) from z to 8r) gives 

Sr] ,dw .d(uw) s , 2r o, n o % /c . \ 

' S— K - — -)dz + £ 6 (nig - u> ) +p s - p= -{5f] - z) 



at 



dx 



Sr, g 



dl 



dx \ dz 



dw 



dw 



+eS / — I Vq— I dz — 2e8vo— — h 2e8vq-— 



0z 



dz 



RR n° 6451 



12 



Sainte-Marie & Bristeau 



and using the boundary conditions (j2"5|) and ((2l?|) . it comes 



£ 2 «5 



d_ 

dt 



w dz + 5 



d_ 

Ox 



5 1 1 



(uw) dz - s 8 w + 5p a - p = -(Sri - z) 



Sr, g 



Ox 



+eS / — — uq— dz — 2e8vo — — 



dw 



Classically we have 



du s 
dx 



du 




^drj du 


du 


dx 


s 


dx dz 


s dx 



dz 



0(e 2 5), 



dz 



(34) 



and using relations (f5T7)) . ((3"4"|) and the Licbniz rule we have 



s5 



Sr, 

dx 



du 

d z 



, „ r dw du c du 

dz - 2ebvQ— = edv — + ebv — 
dz dx dx 



0(e 3 S). 



This leads to the expression for the pressure p 

( d f Sv d [ 

p = Sp a + (Srj - z ) + eH\-J w dz + 5— J 



z 2 8 2 w 2 



w dz + 5— 
dx 

. du du 

sovq- edv — 

ox dx 



(uw) dz j 
+ 0(e 3 S). 



(35) 



Hereafter several models of growing accuracy and complexity will be derived, 
depending on the level of approximation chosen for Eq. (|35[) . In the hydrostatic 
case, we will consider an approximation of p in C(e 2 <5), then in scction[4]we will 
use two expressions of p respectively in 0(e 2 5 2 , e 3 5) and in 0(e 3 8). 



Remark 3 For the derivation of Eq. 133\) note that due to the reseating applied 
to the time derivative of Zb, we have 



5ri du 
. dt' 



dz 



d_ f Sl1 
dt 



u dz 



dr) 

s m Us 



£ dz b 
£0—u b . 
dt 



Remark 4 The second relation in I2ty) is crucial for the derivation of shallow 
water models. When considering large friction coefficients then the assumption 
of asymptotic regime k = ekq no more holds and relation L32\) leads to 

dz 

meaning the assumption of motion by slices has to be justified by other argu- 
ments. 



3.2 Hydrostatic approximation 

We begin with the classical hydrostatic approximation. The objectives of this 
paragraph are twofold. First we want to obtain the expression of u as a function 
of 5, e, Vq, kq and Hg- And second, we aim at verifying that despite the 
parameter 6, we recover the well-known formulation of the viscous Saint- Venant 
system with friction as expressed in the following proposition 
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Proposition 1 The viscous Saint- Venant system defined by 



d{Hu) d(Hu 2 ) g dH 2 



dt 



Ox 



2 dx 



H dp a _ gH ^b_ + d_u uH <^ 



dx dx dx 



1 



3v 



H 



(36) 



(37) 



where H = rj — z b and v = (1, ^§ff) Tu > results from an hydrostatic approximation 
in 0(e 2 S) of the Navier-Stokes equations. 

Proof of prop. Q] wc retain only the terms up to ed in the expression (|35jl for 
the pressure p i.e. we have 



r n ,r \ ? du 8u 

p = dp + (dr)- z) - £ ov — - edv — 



0{e 2 S). 



(38) 



And Eq. {33) with Eqs. {25]) and {38]) gives 



eS 



d(H 5 i 
dt 



d(H s 'i 



s_dHl 
2 dx 



dx 



~e5K u b - e5—(H sP a ) + ed ' g^P a - £ ^Pb + 0(e 2 6), 
that is also using the expression of p obtained in Eq. (1381) 



.d(H s i 



m 



+ <5 : 



,d(H 5 i 



dx 



iafff 

2 dx 



-6K u b - SH S — — pb + 0(eS). 

dx dx 



(39) 



Note that due to the assumption concerning the time derivative of z b and the 
associated rescaling, the first term in the left hand side of (j3"9"]) reads 

d(H s u) du d(rj - z b ) _ 

and (|39|) coupled with (|2"6"|) gives 

.du c o_du T dn Sk c dp a 

5 m +5u dx +5 Tx-~ir s Ub - 5 ^ + 0{£5) - 

Now we come back to Eq. flU, using {2"5"|). (fSS]) and {35]) we get 



S—(u —) 

dz^ dz 



.du 



du dp 



8 



= eS— — h eS u— — h eS w— — h e— e 8— Vq — 



dt 



dx 



dz dx 



di 



dx \ dx 



e8^- + eS 2 u^ + eS^-{r, + P a ) + 0(e 2 S) 
dt dx dx 



u b + 0{e 2 8). 



(40) 
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Integrating from Zb to z and taking into account the boundary condition 
we deduce 



du 
dz 



1 



z - z b 



u b + 0(e 2 ), 



(41) 



and we obtain the following formula which gives an expression of the vertical 
velocity though a parabolic correction 



u = 1 



[z - z b 



(z - z b ) 2 



vo ' 2H S 
Then integrating from z b to St], we obtain 

en 



) u b + 0(e 2 ). 



1 



3^o 



H s u b + 0(e 2 ). 



(42) 



(43) 



Moreover 



which yields 



meaning 



Using (US), (g2J) and 



2£K , (z - Z^ 
Z — Z), 



))^ + 0(e 2 ), 



u 2 = u 2 + 0(e 2 ). 
the right hand side of Eq. 



(44) 



can be written 



£§ d_ f 5ri 
dx 



2v^dz + 6 2 P- P a -^[ Pb + 2eSu ^ 



du 
dx 



-Sk 1 + 



dx 

he 2 ( dzb 
dx 



dx 



du 
dx 



Ub = SKoUb - Hs 



b 

dz b 
dx 



dH s d , Sit, 

p o +e5-5-(2i/ if s -£-, 
ax ax ax 



0(£ 2 <5). 



(45) 



Finally from Eqs. @, J33J), 03]), dHJ) and (gSJ, we obtain the model 
dr] dzb d 



dt dt 



S^^ + S 2 
dt 



dx^ 
d(H s u 2 ) 



dx 



0, 

ldHl 
2 dx 



-e<5— UvqH s — ) 
ax ax 



3i/o ' 

2; 



In terms of the initial variables, the preceding model becomes (|36 |) -(|37 |l that 
complete the proof of prop. [U Note that when the bathymetry is constant 
Zb(x,t) = z®, this formulation is equivalent to the viscous Saint- Venant system 
obtained by Gerbeau et al. Q3] and Ferrari et al. [T3J. 
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4 Two non-hydrostatic shallow water models 

In the previous paragraph we have obtained an approximation of the Navier- 
Stokes equations up to eS terms using an hydrostatic approximation of the 
pressure p. In this section we consider two more acurate approximations of the 
pressure p respectively in 0(e 2 S 2 ) and 0(e 3 5) leading to two non-hydrostatic 
extensions of the Saint- Venant system. 



4.1 First extension, <5 <^ 1 

The first refinement of the classical Saint- Venant model (|36[) - (j3"T)l is achieved 
by considering the pressure p given by Eq. (f3"5|) with the terms up to 0(e 2 5 2 ). 
This means we consider the momentum equation along z is no more reduced to 

dp , d ( . du \ . d ( „ dw\ _ . 2 „ 



dz dx \ dz J ' dz \ dz 

but given by 

and the convective terms are still neglected. Since we keep the terms in e 2 S 
and drop those in e 2 S 2 , this means we assume <5 <gC 1 and due to the applied 
rescaling this implies U -C C so we are in a fluvial regime. The following result 
holds. 

Proposition 2 The system defined by 

f + | («)-«. (46) 

| (ira) + JL (HfiJ) 



d 2 {z b u) 
dxdt 




(47) 



where v = (1, ^r) T 1+ ";if results from an approximation in 0(s 2 S 2 , e 3 5) of the 
Navier-Stokes equations. 

The proof of proposition [5] is given in the next paragraph, we examine here 
some properties of the model (j4^ |) -(|4T |) . 

Note that except for the dissipative terms corresponding to viscosity or fric- 
tion, all the terms added in the non- hydrostatic model (|46|) -(|47 |) compared to 
the original Saint- Venant model (f5u]) - (|3T|) appear as time derivative of the vari- 
ables Zb, rj or u. This means in a stationary regime, the solutions of (|36p ~(137 [) 
and fIB ]) -((I7j l are identical. 
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We first examine the system (pt6l - (|T7|) without friction and viscosity. Starting 
from the Euler equations instead of the Navier-Stokes equations does not allow 
to account for the motion by slices as obtained in relations (|24|) and (|4Tj) . So 
if one wants to neglect the viscosity and friction effects in the model P6]) -P7 ]) . 
it is necessary to consider an asymptotic regime for example under the form 
v = f3v nv , k = f3 2 n n f - and conversely v§ = [3vo,nv> K o = /3 2k o,ti/ - with 
/3 C 1. Introducing the preceding asymptotic regime and considering f3 — ► 0, 
the formulation of {46" ]) -([iT | reads 

+ >* 2 > III 

dz b ( zl d 2 u d 2 (z b u)\ dp a z 2 d 3 z b 

—gti — — „ „ + z b — h - 



dx \ 2 dxdt dxdt J dx 2 dxdt 2 

or equivalcntly in a non-conservative form 

9u _9m c?r/ 9 3 u z b d 3 (z b u) dp a z b d 3 z b 



dt ' 9x ^ ^ dx ^" 6 dx 2 dt 2 dx 2 dt dx ' 2 dxdt 2,1 

that is analogous to the expression obtained by Peregrine [20] ■ It is worth being 
noticed that, in any case, the formulations obtained by Nwogu [T5], Walkley [23] . 
Saut ei aL [5] and Soares Frazao et a/. |22j are different from the preceding 
ones. The differences lie either in the continuity equation or in the momentum 
equation. 

The mathematical and numerical analysis of the obtained model is not in the 
scope of this paper but let us mention some interesting works in the literature. 
The Sobolev equation 

d d 2 u du d . . .du du 

~dx- {a{x) dxd-t ] + c{x) m = -di^dJ + m dx-> 

has been studied by several authors [U 0] as an alternative to the Korteweg- 
de Vries equations. Perotto and Saleri [21] proposed an a posteriori error 
analysis for the Peregrine formulation of the Boussinesq system with constant 
bathymetry. Bona et al. [5] [S] have studied the well-posedness of several high- 
order generalizations of the Boussinesq equations. 



4.2 Derivation 

Proof of prop. 0' the refinement of the classical Saint- Venant model ([SBf - lpT?]) 
is achieved by improving the approximation for the pressure p. Actually, if wc 
only drop the terms in 0(e 2 6 2 ) in the momentum equation along z so the system 
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- (fTU|) becomes 

dzb d 
dt dx 



u dz, 



c <9m ^ 2 du 2 2 

^ + £ = -1 + 

at c»z 



duw 


dp 


e2 4 


dz 


dx 


dx 


~( 

dx V 


e5vq 1Tz) 


dz 



du 

dx J <9z 
az 



9z 



with the boundary conditions (|28[) - (|32p . This means we consider the pressure p 
is given by (f3l>|) where we retain only the terms up to e 2 ^ 2 and e 3 S i.e. 



5p a + (5r) - z) - s5v a 
+0(e 2 8 2 ,e z 8), 



da 
dx 



eSv 



du 
dx 



.d f Sr > 



w dz 



leading to 



Pnh = 8p a + (Srj - z) - eSv ^ 

dx 



du d 2 Zb 

e6 „ _ +£S (5ri-z)-^ 



d t ' d r z 
5 lFt J ikj udz i dz + °( £ 6 > e s )- 



(48) 



Retaining only the terms up to 0(e 2 5 2 , e 3 S), relation ([55)) gives 



S dt 



d r Sr > 



d r Sn 



u dz + 8 — / u dz + — 



d r Sl1 



dx 



dx 



611 . du 



dx 



Pnh dz 



dr) a dz b 



s8— / 2v Q —dz + 5 z -±p a + — - 2e8v Q — 



dx 



dx 



dx 



dw 



dz 



Pnh\ b 



(49) 



Now we derive the expressions for the quantities appearing in (|48[) and (|49[) and 
depending on u, w and p. Since k — « 0i i + 0(e), from Eqs. (jUJ) and (|4*3"|) wc 
have 



du 

dx 



ekq , (z- z b ) 2 A du b 
1 H [z — z b 1 



9a; 



£Ko / <9zf, 



\ dx 



-1 



Z - Zfcx 9ff«5 (z - Z 6 )' 



1 i £K ° ( 

1 H (2 — 2& - 

^0 



#5 ' 

(z - Zf,) 2 \ du b 



u b + 0(e 2 ) 



2H S 



dx 



ekq dz b 
vq dx 



z - z b (z - z b f 



Hs 



2H 2 



u b + 0(e8), 
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so 



du 
dx 



du 
dx 



6fo ' dx Qvq dx 



= (1 

du b ekq dzb 
dx vq dx 



u b + 0{e 2 ) 



(50) 
(51) 
(52) 



'l £K ° U\ 9U £K ° ( dZb , ldH S\- , ml 2^ 

^-^d-x-W { ^ + ^ )U + 0{£ } ' ^ 



and 



8q 



du , EK 

vo-z- = v Hs[l + - — 

dx 6Uq 



H s du b ek H s f dzb 1 dH s 



dx 2 v dx 3 9a: 

du ekq dHs _ £kqHs ,dz b 1 dHs 



d r 



3 5 dx 



= U ° Hs dx 2 
and finally from (|34[) we get 



du _ £ ^i Hs ( ^1 + 1 ^E£ \ 



dx 3 dx J 



dx 6 dx 
u + 0(e 2 ), 



)u b + 0(e 2 ) 
)u + 0(e 2 ) 



(54) 



Sri du 
dx 



6,1 du 



dx 



EK 



= v Q Hs 1 + —1H S — 



6v 



du skq dHs 



dx 



6 dx 



H s u + 0{e 2 ). 



From ([51]) and |53|) we have 



Ph \ b -2eSu — 



5p a + H 5 + e5v |^ 
dx 



. du 
ebvQ — 
dx 



2 K TT du 



= Sp a + H 5 -e z S^-H s 

2 dx 



-s 2 5k 



dz b | ldH s \ 
dx 2 dx J 



u + 0{e 3 S), (55) 



and 



5r] 



d 

(2e^ ^ ~ Ph)dz = ~H s 5p a - ^± 4 



2 

Hs du 7 dz b _ 6 drj 



du 
dx 



6 9i 6 c?o; 



3 9x 



(56) 



where corresponds to the gravitational, viscous and friction part of the pres- 
sure p given by Eq. |4*5|) i.e. 



Ph = 5p a + (6r) - z) - eSv ^ 

dx 



ax 



(57) 
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Inserting (|55[) and ([5(J| in equilibrium ((19"]) leads to 

-H s —(Sp a + z b ) + e5—(^H s —) 
e 2 S d ( ( du _dzb _\\ e 2 5n dzb ( dzb _ du 



6 m \ \ ox dx J J 2 dx \ dx ox 



where Ap n h = p n h — Ph- And using the expression for the pressure p n h given in 
Eq. (0E1) it comes 



i Apnhdz = - £5 ^ + Zb) dxm +£6 ^t^bT 



( Sr,^ - °<§®) + e 2 6^-^ (58) 



dt \ 'dx dx J 2 dt 2 

, z 2 ( z b d 2 u d 2 (z b u) d 2 z b \ ,~, 2X 2 _ 



and 



We finally obtain the model 

^ + A (ffifi)=0 , (60) 
^ i tt —\ . r2 9 m _ 2 , , i 2 a / z b 3 9 2 s z 2 a 2 (z fcM ) 



<5-(i^) + ^(W) + - — A. - s 2 S— 



dt dx 2 dx dx \ 6 dxdt 2 dxdt 

tt ° ,c a ^ ° ( a e tt SKq ( du dZb _ 

-H s — (5p a + z b ) + — 4e<5j/ ffi- T - + -r^b 2 b — + 7— u 
ox dx \ dx 6 \ ox ax 

e 2 Sk dzb ( dzb _ du\ 9c dzb ( zi d 2 u d 2 (zbu) 
u + z b — + e o—- --f — — + z & - 



2 dx \ dx dx J dx \ 2 dxdt dxdt 

-** (i + * - A IS + < 61 > 

that complete the proof of proposition^ When the terms in 0(e 2 S) are dropped 
in (l(il[) . we verify that we recover the classical viscous hydrostatic Saint- Venant 
model with friction (|55|) - (|57|l . 

4.3 Energy equality 

Until now, we have not verified the derived models satisfy an energy equality. 
The system (p^ -lfST ]) that is equivalent to the Saint- Venant system, admits a 
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dissipation energy [3J [7] . Indeed we have 

dE h d ( , H 2 . du\ dp a ,du. 2 

I u(E h + g—) - 4vHu- = -H^- AvH{-) 



dt dx \ v 2 ' dx J dt K dx 

k(v,H) - 2 . rjdz b 

with E h = Sf- + a H ^+ z ") + Hp a . The energy equality ||£SJ) associated with the 
hydrostatic Saint- Venant model can be obtained using classical computations 
by multiplying Eq. (|33|) when p = ph by the velocity u. 

The only differences between the hydrostatic Saint- Venant model ([3T>)) - ([3"T[) 
and its extended version PB| - (jT7)) comes from 

• the non hydrostatic terms of the pressure p n h, 

• the terms involving the viscosity and the friction at the bottom, 

so the energy equality for ([50)) - ([CT|) will differ from Eq. ([55)) only by the terms 
_ d f Sv _dz b 

Cl = U— J Ap nh + U— Ap nh \ b , 

d f Sri (_ . du 



C2 = n Tx L 2e ^>' 



d_ f Sr > 
" dx . 



C 3 = U— I p vJ , 



r -dz b ( dw 

Ci = u— 2edv — 
ox \ dz 



dz b du 



where Ap n h = p n h — Ph and p v j = p b — 5p a denotes the terms in the pressure p 
containing the viscosity and friction. The quantities C1-C4 corresponding to the 
non-hydrostatic terms, come from the multiplication of Eq. (|33[) by u and have 
to be added to $2$. Since u = u + 0{e) = u b + O(e) and Ap nh = 0(s 2 S 2 ), we 
rewrite C\ under the form 

d f Sv dz h 
Cl = U dlc J ^Pnh + u b — A Pnh \ b + 0{e 2 8 2 ) 

8 f 5ri f 5ri 8u A dz b /n . 2r2 , 

= — J upnh- J T^Ap nh + u b — Ap nh \ b + 0(e d ) 

= A p uApnh + [ W Ap nh ] 5 Jl - p + u b ^ A Pnh \ b + 0(s 2 S 2 ) 

d f 5r > A A [ 5r > 8A Pnh dz b 



dx 



J " uA Pnh + w s A Pnh \ s - J V w 9 -^- - ^ A Pnh \ b + 0(s 2 S 2 ), 



1 z b J z h 

where relation ([3"Tj) has been used. From Eqs. and ([Ff]) . we have 



dt 



Ap nh \ s 


= 0(e 2 S 2 ), 




r Sr/ 


APnh\ b 


= e 2 6 / 




J z b 


dAp nh 


2 ?dw 


dz 


= - £5 m 
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leading to 

Due to the rescaling applied to the time derivative of z b (see paragraph 12. 2|l . 
the Leibniz rule applied to obtain the preceding relation reads 

f 5ri dw ,dz h wl ,dr]w 2 s d f Sr 'w 2 d f 5ri w 2 

And finally we have for C± 

d f Sv a 2,9 [ Sr > w 2 8z b m . 2 ^ 

Ci = ^ uA Pnh + e 2 S-j Zb ---A Pnh \ b + 0(eW). 

From relations and we obtain 
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so we have for and C3 



d 



dx 




Si) 



-I «l»o S 



du 
dx 

du 



S V ° L. b dx \dx 



d_ 

dx 

d_ 

dx 

d_ 

dx 



zb 

8r) 



Up v j 



dx 

5t > du 

, dx 



6V fHHY E ^L f 6V d ( fU ^ 
dx / dx /, t dx 



0(e 3 S), 



Pv,f, 



Sri dw 



S 77 

Up v + [Wpvj] 5 ^ - sSuq 

b 

5 " d(fu) 



Si) 



Si) 



d(fu) 
dx 

d 2 u 
dxdz 



dx 



'- PvJ + 0(e 6 S), 



d f Sn _ 

fa J UP V J + [W PV J] 



S V Q 



PvJ + 0(e 3 S), 
— 2e8v 



dw 



dz 2 



2e8v / — w— + 2e8v / — 



dz V dz 



5r > (dw 



-eS i/ 



dx 



5 rj 



du 
dz 



c2 drj du 
eb v —w s — 
dx dz 



c dz b du 
edvQ—w b — 
dx dz 



. [ 5r > dw du 
Jz b ox dz 

d [ Sr > _ 
fa J up v j + [wp v j] 



(>// 



Zb 
(>// 



d(fu 
dx 



PvJ +0(e 3 8), 



Sv 'dw 
dz 



- e8v 



du 



dw 



dx 



61 j 



dz 



-2 drj du 
ed vo-^-w s — 
dx dz 



. dz b du 
-Edi> —w b — 
dx dz 



£ V ° Jz b dx dz 



Si) 



d(fu) 
dx 



PvJ + 0(e 3 8), 



and from relation fill ) we also have 

1 



z b 



Sr > dw du 
dx dz 



Zb 



5r > dw 
dx 



z - z b 



u b + 0(e 2 ), 



EK - 



H s d 2 z b 
2 dxdt 



U b + SKq 



St} 



d 2 u d 2 {z b u) \ 
dx 2 dx 2 J 
H s d 2 z b H 2 d 2 u H s d ( dz b 

£K0 ^dxm Ub - e ^dx- 2n + £ ^ 



2 dx \ dx 



z - z b 
+ 0(e 2 



u + 0( 



H s d 2 z b 

£K °—dxm Ub - £Ko dx- 



d (H}du 
— i 

6 dx 



en 



Hi 



du 
dx 



Hs dHs du _ 

+EKq — — — — U + EKq - 



d f H s dz b _ 2 \ en dHs dz b _ 2 

r.n II. u. 



3 dx dx ' dx \ 2 dx 
The preceding expression shows that due to relation (jUJ), the term 

" Sn dw du 



2 dx dx 



0(e 2 



dx dz '' 
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has to be treated as a friction term in the energy equality. We finally have for 
K = C 2 -C 3 + C 4 

d [ Sr > / r du \ n c f Sr > ( fdu\ 2 fdw\ 2 " 

+s Su— (w—\ + e 2 6— (^—u-^u 2 

dx J Zb \ dz J dx\ 2 \ 3 dx dx 

e 2 Sn a ( + 9H S _\ 2 _ /_ 2 / dzb\ 2 + s dr) dz b 

6 \ V dx dx J \ \ dx J dx dx 



Returning to the initial variables and integrating C\ and 1Z into relation 
gives an energy equality for the model (fi^ |) - (|T7|) under the form 

d / Hvp\ d ( ,„ TT . fV„ TT <9u <9u 

a + — J + to r {Eh + Hpnh) - " 1 1 2 ^ ^ + ^ 

+ — (k ( ^-—u + ——u 2 

dx \ \ 6 to 2 9i 



~ 2vl 11 toj + toJ rehto 1 to 



11 /to, x 2 



dp a ( du b \ dz b z b d 2 z b _ 

-H— + Pnh \ b + 2v— \ — + K-^-u, 



where 



Hw 2 = w 2 



n ( dz b du | d(z b u)\ z 2 f du 



v dt dx dx J 3 \ to 

/ d(z b u)\ 2 2 dud(z b u) ( dzh~ 2 

- Zb [-^r) +Zb d-x—~ Zb \-dt 

+ b ~dt \ 2dx~~~dx~ 

Hp nh = / p nh . 

J z b 

When the time derivatives of p a and z b are dropped, the right hand side of the 
preceding energy equality is always negative. 

4.4 A more complex approximation, 5 = 0(1) 

Now we return to the dimensionless and rescaled variables. The assumption that 
the elevation of the free surface is small done in paragraph 14.21 is now relaxed 
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i.e. S = 0(1). This means that no assumption is made concerning the hydraulic 
regime. We consider for the pressure p the complete expression obtained in (j3"5]) 
and the following proposition is a refinement of the Proposition [2] 

Proposition 3 The system defined by 

— (Hu) + —(H u 2 ) + + 9 ( g P"g-"») = -H— - gH— 

dt dx m 2 dx dx dx dx 

+ ® (AvH®^) + ® ( k(v H)h( H ^ U (ldz b ^\dr\ 
fa ^ dx' dx \ ' \ 6 dx \ 6 dx 3 fa 

n(v b ,H) ( 5 (dz b \ 2 \ _ dz b f A dH dz b ._ 



1 I "(vi-.g)g I 2 V dx J / v ' ; fa V v 2 dx dx 
H du\ dz b dz b d 2 z b Id ( 2 d 2 z b \ 

+ T fa J - & p " I* + Zb fa ^ " 2 fa ^ a^J ' (64) 

where v = (1, ^r) T results from an approximation in 0(e 3 5) of the 

Navier-Stokes equations. In the previous expressions, H m is a modified wa- 
ter height taking into account the Coriolis-Boussinesq coefficient and p ng ,nv, 
Png,nv\ b corresponds to the vertically averaged and bottom value of the non grav- 
itational and non viscous part of the pressure p given by 135\) . 



Proof of prop. [H' we still start from the averaged momentum equation (|33[) 
where, compared to the first extension of the Saint- Venant model detailed in 
paragraphs 14.11 14.21 and 14.31 the expressions of 

p dz and / u 2 dz, 

z b J Zb 

have to be refined. The approximation u 2 = u + 0(e 2 ) obtained m para- 
graph is no more sufficient. From d3SJ), (@2J| and (g3J) we get 

H s ))u + 0(e 2 ), 

w = ^-i-[[ z - z b + £ f(z,z b ,H s )dz)u) +C(e 2 ) 





e«o 


Vq 




ef(z 


dz b 


d 


dt 


dx 


dz b 


d 


dt 


dx 






5— 

dx 


(P a + 



Zb 



2\ 



;, z b , H s )u) + (D(e 

and Eq. (j4l)|) is equivalent to 

c du C 2-du -fa J „ Sk 

S^T + S 2 U— + S— 1 - = ~S—p a -. r 7i 

dt dx dx dx 1 H s (l + f^H s ) 

+ ^fa^fa) + ° (£ ^- 
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Now we can improve the approximation (|40|) in the following way 

. d , du-. T du du r9 du dp , . d ( du\ 

dz dz at dx dz dx dx \ ox J 



+ ^<" + "°) + ^£- 3rt s(™£) +o < rt ) 



H s (l + ^H S ) dz dt yj ' dx 

_ 2 () / du\ 2 4u dHs du 3 



^ c?a; \ 9xi ^ £ ^ Hg dx dx^^ E 



Taking into acount the boundary condition (|32|) . an integration of the preceding 
relation from Zb to z gives 

9u £Kq N 2 — Zb\ 2 sr f c- ji ?-du 



= (iTj^H 1 --fjrJ 5+e H {a/ '" + ^X. / 

♦iK') + '£(*jC') + <*-*>£(*£ 

where the relation 

J w— = eu J w— = eu{fw - J\bw\ b ) + eu— J ,/, 

has been used. Another integration of relation (|65p between Zb and z gives 

1 H 2-26 H u / fw 

vo V 2H S J J v J Zb 

e 2 6 2 _du ( z f zi „ e 2 S d (_ f z f Zl „\ e 2 S 2 d ' rZ rZl 



dx J Zb J Zb v dt V J Zb J Zb J vq dx 
e 2 S , , 2 d 2 u 2 .2(z - z b ) 2 dH s du _- 3r . 

l + * + ■•**. + <><.■*>. 

so we obtain the new expressions for u, u 2 and u 2 



1 + -^H s u b + e'SAu + 0{e A 8), 
_ I i i 2eK ° tr. \ „.2 



/ + — irJul I f)+—i-lu 2 l I f 



u z = \l + -^-H s jui + 2e 2 SAu + 0{e 6 6), 
2 / 2en Q ( (z~z b ) 2 \ e 2 n 2 ( (z - z b ) 2 

u = [ 1+ — { Z - Zb — 2iir) + — { z - Zb -^r 

+2e 2 6Au + 0(s 3 5), 
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so finally 

"5 — Hs + TH - 



1 + ^#,5 + t ^fH 2 s ) u 2 + 2e 2 (5Au 



1 + &H 2 ) u 2 



0{e 3 S). 



Now concerning the expression of the pressure terms, it has to be noticed that 
Eqs. (f3l>|) and l[48p only differ by the terms 



Using 



d r Sv 

A = e 2 S 2 ^- I uwdz-e 2 6 2 - 
dx 



. . dz b d 

u = u + 0{e), w =-Qj--g^l udz ' 



Zb 



it comes 

_2r2 / u y - z u , TO, ^ 2 



£ 2j2 / S 2 rj 2 — z 2 d ,_du, ^ 2 <9?7 cht_ ^ d ,_d{z b u) , 

' dx^ dx' dx dx dx ^ dx • 



c dr] d(z b u) _ ( du d(z b u) 
o— u — —z 



dx dx \ dx dx 

This leads to the new expression for the fluid pressure p appearing in (|33p 

p dz = / (p nh + A)dz 



zi, J z b 



5ri , e 2 5 2 H s ( ATT2/ du. 2 nTT2 d 2 u nTT dH s du 
Pnhdz + ^ ( -4if 5 2 (— ) - 2Hju— - 6H 5 - 



6 \ v dx' dx 2 dx dx 

dz b du d 2 z b 9 dzbdHs 9 \ m , q„ n 

+§H s -±—u + 'AH S — ±u 2 + 6-77— —p——u + 0(e 3 S), 
dx dx dx 2 - dx dx 



where i v p n h = L. (ph + ^Pnh) is given by (J35J). Conversely using we 
obtain 



where p T i/i| b is given by Inserting ([5l|) . |55|) and (j56"|) in equilibrium (j3"3"|) 

leads to the system 

, I (JM) + s >jl IHi ,^ } + \m + . + 



^^(4^ # ^H) + e 2 jiL f w Hs( — — -(-— + 



u 



2r <9z b / AdH s dz b ._ H s du\ dz b 
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where 



2^2 



HsPng,nv = (P ~ Ph) dz 

J Zh 



d /iff d 2 u H 2 d 2 {z b u) H$ d 2 u 



£ 5 dx\ 6 dxdt + 2 dxdt Sv 2 dxdt 

-<(^-^))^(--«£) 2 

„„ 9 <9 2 u „„ dHs du „ TT dzb du 
ox 1 ox ox ox Ox 



i9x 2 dx dx J 2 9i 



and 

Pn fl ,m>| 6 = (P-Ph)| 6 



dt 2 dt dt J dx dx dx dx 



dxdt' J \ dt 2 dt dt J dt 2 

e 2 8 2 ( d , TT ndu _.\ 9 , dz;, 



In terms of the initial variables, the model (f66l)- (167D corresponds to the one 
depicted in proposition[3]with obvious expressions for H m , Hp ng ^ nv and p n g,nv\ p - 
In order to obtain the energy equality for the model (f63|) - ([64|) . we use the 
same process and the same notations as in paragraph s. 31 but the approximation 
order is now 0(e 3 S) instead of 0(e 2 S 2 ). Still using u — u + 0(e) = u b + 0(e), 




-dz b 



[ Sr > dAp dz b , nl 



with Ap = p—pnh, P being given by |55|) . From Eqs. (|3"5|) . (|5?|) and the boundary 
condition (|28j). we get 

Ap| s = 0( e 3 <5), Ap| b = £ 2 <5 2 1 5 " ^ + £ 2 5> 2 - w 2 ) + 0{eH), 
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leading to 



dAp ?T? d(uw) „ ,., dw 
— - = -e 2 5 2 ^—^- - 2e 2 d 2 w— + 0(e 3 6), 
oz ox oz 

d f 5 " . 0.0 f 5n duw 2 2r2 , , 

Zb J Zb 



C X = - uA P + e^ W — + -s^( W i- wt ) 



Ap\ b + 0(e 3 6), 



dz b 
dt 

d ( f Sv /A , w 2 ^\ dz b 
dx 



J\ { A P + ^)\-^A P \ b + 0(e 3 S). 



Returning to the initial variables, the preceding relation and the expression of 
1Z obtained in paragraph 14.31 allows us to write an energy equality for the model 
(ED-lES under the form 



-{E + Hp)- V (2Hu-+ W -)+ K (—-u-^^u 2 



dt dx V 7., V dx dz J \ 6 dx 2 dx 

. du\ 2 / du\ 2 \ k ( jj9u dH _ 
dx J \dx J J 6 \ dx dx 

k I ( dz b 1 dri\ 2 If dri\ 2 \ _ 2 k ( 3 / dz b \ 2 \ _ 7 

u - -w 1 + - hr- w 



with 
E 



Z\\dx AdxJ 8\dxJ / 1 + #M 2 V 9a: 

<9p° / <9u b \ 9z fc H d 2 z b _ 

-H— + l Pnh \ b + 2v— \ — k-—u, 



Hu 2 Hw 2 gH(rj + z b ) 



2 2 2 



15i/ 2 



^ = [\> = H( r]2+ ^ + f g V - + z b ) dU dM 



3 \dx J dx dx 

d{z b u)\ 2 \ i jj( 9zb \ 2 dZb ( r,2 ~ z b dU i H d ( z bU) 



dx J I \ dt J dt \ 2 dx dx 

Note that except for the friction terms, the previous expression is analogous 
to the energy equality for the Navier-Stokes system [16] but expressed with 
the vertically averaged variables. When the time derivatives of p a and z b are 
dropped, the right hand side of the preceding energy equality is negative when 
^ is enough small. 



5 Conclusion 

In this paper we have derived two extensions of the Saint- Venant system when 
the hydrostatic assumption is relaxed. The obtained models, especially in sec- 
tion [H are similar to Boussinesq type models but derived in a more rigourous 
context and satisfying an energy equality. 
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On one hand the averaged models of shallow water type presented in this 
paper reduce the complexity of the discretization of the Navier-Stokes equations 
since they are written over a fixed domain. But on the other hand their math- 
ematical formulation is more complex since high order derivatives - especially 
in space - appear. 

The preliminary numerical simulations and comparison with experimental 
measurements performed with the proposed models are promising. They are 
not presented in this paper and will be described in a forthcoming publication. 
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